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Using Ha as a Tracer of the Emission Region of LS I +61 303 

M. Virginia McSwain 
Lehigh Umverstty, Bethlehem, PA 18015, USA 

The 7-ray binary LS I +61 303 is one of the brightest Fermi sources, with orbitally modulated 
emission across the electromagnetic spectrum. Here we present Ha spectra of LS I +61 303 that 
exhibit a dramatic emission burst shortly before apastron, observed as a redshifted shoulder in the 
line profile. A correlated burst in radio, X-ray, and GeV emission is observed at the same orbital 
phase. We interpret the source of the emission as a compact pulsar wind nebula that forms when 
a tidal mass stream from the Be circumstellar disk interacts with the relativistic pulsar wind. The 
Hq emission offers an important probe of the high energy emission morphology in this system. 
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I. INTRODUCTION 

LS I +61 303 is a high mass X-ray binary (HMXB) 
that consists of an optical star with spectral type BO 
Ve and an unknown compact companion in a highly 
eccentric, 26.5 day orbit [3 . While the system has a 
relatively low X-ray luminosity for a HMXB, LS I +61 
303 is the 15th brightest 7-ray source included in the 
Fermi LAT 1-year Point Source Catalogue ([5]). The 
Be disk interacts with the compact companion, pro- 
ducing emission that has been observed to vary with 
orbital phase at every wavelength across the electro- 
magnetic spectrum, from radio to TeV (eg. [T], [101 )• 
[lOj found periodic radio outbursts that peak near 
(/)(TG) = 0.6 — 0.8, and they defined the arbitrary 
reference for zero phase at HJD 2,443,366.775 that re- 
mains the conventional definition for LS I +61 303. 
Periastron occurs at (/'(TG) = 0.275 j3j. 



II. OBSERVATIONS 

During 2008 October and November, we performed 
an intense multiwavelength observing campaign on LS 
I +61 303 supported by a Fermi Cycle 1 program. 
We obtained optical Ha spectra of LS I +61 303 at 
the KPNO Coude Feed telescope over 35 consecutive 
nights to study the evolution of the emission during a 
complete orbit |H]- The Ha line profile exhibits a 
dramatic emission burst near 0(TG) '-^ 0.6, observed 
as a redshifted shoulder in the line profile (see Fig. 
[1]) as the compact source moves almost directly away 
from the observer. 

Smaller temporal changes in the red spectra sug- 
gest additional Ha emission variability, so we sub- 
tracted the mean emission line profile to investigate 
the residuals carefully (see Fig. [2]) . During about half 
of the orbit, 0.9 < (/)(TG) < 0.6, the difference spectra 
reveal a partial S-shaped pattern similar to a spiral 
density wave that is commonly observed in Be star 
disks [S]. [H] also observed a strong blue peak near 
(/>(TG) = 0.23, which supports the development of a 
spiral density wave near periastron. After this phase, 
the peculiar red shoulder develops. 

We measured the equivalent width of Ha, M^hqj for 



each spectrum by directly integrating over the line 
profile. (We use the convention that Whq is negative 
for an emission line.) The errors in Who are typically 
about 10% due to noise and placement of the contin- 
uum. Figure [3] shows that during our Coude Feed run, 
W]ia decreased slightly just before periastron. Since 
Who is correlated to the radius of a Be star's circum- 
stellar disk '5], we interpret the decline in emission 
as a slight decrease in disk radius as gas is stripped 
away by the compact companion. Whq then rises dra- 
matically with the onset of the red shoulder emission 
component near (/)(TG) ~ 0.6. 

Figure [3] also compares our recent Who with those 
measured by [3] . Their data were accumulated over six 
different observing runs over 1998-2000, and the long 
term differences in emission strength are substantial. 

Also during 2008 October and November, G. Poo- 
ley obtained nearly simultaneous radio flux coverage 
with the Arcminute Microkelvin Imager (AMI) array. 
The 15 GHz AMI light curve (Fig.|4]) reveals emission 
that peaks at the same time as the Ha "red shoulder" 
outburst. Contemporaneous RXTE light curves from 
[Tljand Fermi light curves ([4]) also reveal orbitally 
modulated emission that peaks just before the Ha red 
shoulder, although their wide phase bins may mask a 
true correlation. The Ha emission clearly traces the 
high energy emission region in this system. 



III. CONCLUSIONS 

The unusual broadness of the Ha red shoulder emis- 
sion is consistent with a Balmer-dominated shock 
(BDS; [7j)- BDS are traditionally observed around 
supernova remnants but are also sometimes produced 
within pulsar wind nebulae and other evolved stellar 
systems. They form when high velocity (200-9000 km 
s~^) shocks collide with the interstellar medium, man- 
ifesting themselves as optically emitting filaments. 
Energetic particles and/or photons may be gener- 
ated in the post-shock region of the coUisionless, non- 
radiative shock. Direct coUisional excitation of the 
pre-shock atoms produces a narrow emission line com- 
ponent that reflects thermal conditions within the pre- 
shock gas. If the energetic particles exceed the shock 
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FIG. 1: The upper plot shows the Ha hne profile over 35 continuous nights of observation at the KPNO Coude Feed 
telescope. The lower plot shows a gray-scale image of the same line. Note that the lower plot of gray-scale spectra are 
not folded by orbital phase but are placed in the same chronological order. From [S|. 



velocity, the pre-shock hydrogen atoms also exchange 
electrons with post-shock protons, manifesting them- 
selves as broad neutral hydrogen lines (widths ^ 1000 
km s~^). The Ha line structure in LS I -1-61 303 is 
complicated by the superposition of emission from the 
circumstellar disk; however, the broad red shoulder is 
consistent with such a BDS. The temporary nature of 
the red shoulder, as well as the correlated GeV-radio 
emission, suggests that the BDS only forms when a 
high density tidal mass stream interacts with a pulsar 
wind in LS I +61 303. 
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FIG. 2: The Hq difference spectra of LS I +61 303, shown in the same format as Figure 1. Note the distinctive peak in 
emission that appears to follow a partial S-curve in the spectra. A fast moving "red shoulder" of emission produces a 
strong redshifted feature between (^(TG) = 0.6 — 0.7. From [S]. 
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FIG. 3: The emission strength of Ha over muhiple orbits is indicated by the equivalent width of the hne (Wua). 
Measurements from the six runs of Grundstrom et aL (2007; identified by the date HJD-2,450,000) and the mean Wna 
from each night of our 2009 WIRO run are shown folded by orbital phase. Note the consistent peak near 0(TG) ~ 0.6. 
From [H]. 
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FIG. 4: 15 GHz radio flux of LS I +61 303 from the AMI array during the dates corresponding to the 2008 Ha 
spectroscopy campaign. Data courtesy of Guy Pooley. 
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